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_J Summary. We study a Life-like cellular automaton rule 52/5*2345 where a cell in 

. £3 state '0' takes state T if it has exactly two neighbors in state 'F and the cell remains 

( — i in the state 'F if it has between two and five neighbors in state 'F' This automaton 

is a discrete analog spatially extended chemical media, combining both properties of 
sub-excitable and precipitating chemical media. When started from random initial 
configuration B2/S'2345 automaton exhibits chaotic behavior. Configurations with 
low density of state 'F show emergence of localized propagating patterns and sta- 
tionary localizations. We construct basic logical gates and elementary arithmetical 
circuits by simulating logical signals with mobile localizations reaction propagat- 
ing geometrically restricted by stationary non-destructible localizations. Values of 
Boolean variables are encoded into two types of patterns — symmetric (False) and 
asymmetric (True) patterns - which compete for the 'empty' space when propagate 
in the channels. Implementations of logical gates and binary adders are illustrated 
explicitly. 



1 Introduction 

Computational universality of Conway's Game of Life (GoL) cellular automa- 
ton [14] has been already demonstrated by various implementations. Most 
famous designs are realization of functionally complete set of logical func- 
tions [40], register machine [7], direct simulation of Turing machine [9, 39], and 
design of a universal constructor [18]. These implementations use principles 
of collision-based computing [7, 2] where information is transferred by local- 
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izations (gliders) propagating in an architecture-less, or 'free,' space. The the- 
oretical results regarding GoL universality is only the first half-step in a long 
journey towards real-world implementation of the collision-based computers 
as unconventional computing [42, 1]. Controllability of signals is the first ob- 
stacle to overcome. Despite their mind-whirling elegance and complexity-wise 
efficiency of implementation the 'free-space' computing circuits are difficult 
to fabricate in physical or chemical materials [3] because propagating local- 
izations (solitons, breathers, kinks, wave-fragments) are notoriously difficult 
to manipulate, maintain and navigate. 

The easiest way to control patterns propagating in a non-linear medium 
circuits is to constrain them geometrically. Constraining the media geomet- 
rically is a common technique used when designing computational schemes 
in spatially extended non-linear media. For example 'strips' or 'channels' are 
constructed within the medium (e.g. excitable medium) and connected to- 
gether, typically using arrangements such as T-junctions. Fronts of propagat- 
ing phase (excitation) or diffusive waves represent signals, or values of logical 
variables. When fronts interact at the junctions some fronts annihilate or new 
fronts emerge. The propagations in the output channels represent results of 
the computation. 

The geometrical-constraining approach is far from being graceful (par- 
ticularly, comparing to collision-based paradigm [2, 23, 34]) but simple and 
practical enough to be used by experimental scientists and engineers. The ge- 
ometrical constraining is successfully applied in design of several laboratory 
prototypes of chemical (all based on Belousov-Zhabotinsky reaction) comput- 
ing devices: logical gates [43, 41], diodes [25, 13, 36], counters [22], coincidence 
detectors [15], and memory [36]. 

What members of GoL family offer us most 'realistic' approximation of 
dynamical processes in spatially extended chemical computers? In its origi- 
nal form, GoL automaton is a discrete analog of sub-excitable chemical me- 
dia [8, 24]. Localized wave-fragments (reaction-diffusion dissipative solitons) 
in sub-excitable Belousov-Zhabotinsky reaction [8, 24] are represented by glid- 
ers [14, 45, 16, 33, 38, 4] in GoL cellular automaton. See examples of direct 
comparisons between experimental chemical laboratory results and cellular 
automaton models in [12]. 

There is also a family of Life-life rules, where 'cells never die,' the state 
'1' is an absorbing state. This is the family of Life without Death (LwD), 
invented by Griffeath and Moore in [16]. In the LwD automaton we observe 
propagating patterns, formed due to rule-based restrictions on propagation 
similar to that in sub-excitable chemical media and plasmodium of Physarym 
polycephalum [5], but no complicated periodic structures or global chaotic 
behavior occurs. The LwD family of cell-state transition rules is an automaton 
equivalent of the precipitating chemical systems. This is demonstrated in our 
phcnomcnological studies of semi-totalistic and precipitating CA [4], where 
we selected a set of rules Life 2c22, identified by periodic structures [35]. The 
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clans closest to the family 2c22 are Diffusion Rule (Life rule B2/S7) [27], all 
they also into of a big cluster named as Life cfc22. 5 

In present chapter we study a Life-like cellular automaton, which some- 
what imitates properties of both excitable and precipitating reaction-diffusion 
chemical systems, and who how to implement a sensible computation in such 
type of cellular automata. We develop a model where the channels are con- 
structed using static patterns and computation is implemented by propagat- 
ing patterns of precipitation. The results are based on our previous studies of 
reaction-diffusion analogs 52/523456 [26] and B2/S2345678 [28]. 

2 Life rule B2/S2345 

Life rule B2/52345 is described as follows. Each cell takes two states '0' 
('dead') and '1' ('alive'), and updates its state depending on its eight closest 
neighbors: 

1. Birth: a central cell in state at time step t takes state 1 at time step 
£ + 1 if it has exactly two neighbors in state 1. 

2. Survival: a central cell in state 1 at time t remains in the state 1 at time 
t + 1 if it has two, three, four or five live neighbors. 

3. Death: all other local situations. 

Once a resting lattice is perturbed (few cells are assigned live states), 
patterns of states 1 emerge, grow and propagate on the lattice quickly. 

A general behavior of rule -B2/S , 2345 can be well described by a mean field 
polynomial and its fixed points. Mean field theory is a proved technique for 
discovering statistical properties of CA without analyzing evolution spaces of 
individual rules [21, 30, 10]. The method assumes that elements of the set of 
states are independent, uncorrelated between each other in the rule's evolution 
space. Therefore we can study probabilities of states in neighborhood in terms 
of probability of a single state (the state in which the neighborhood evolves), 
thus probability of a neighborhood is the product of the probabilities of each 
cell in the neighborhood. 

The mean field polynomial for rule .62/52345 is as follows: 

Pt+i - 28 P Ut + ^pUt + WpUt + ™pUt + WPtll (i) 

High densities of regions dominated by state 1 correspond to maximum 
point near p — 0.35. The average density is reached with one stable fixed point 
p = 0.46 when automaton find its stability around of 37% of cells in state '1.' 
Some interesting behavior can be found in extreme unstable fixed point when 
p = 0.04 (complex class, see mean field classification in [30]). Looking on 
configurations with less than 4% of cells in state '1' we can observer gliders, 

5 http : / / uncomp . uwe . ac . uk/ genar o/Lif e_dc22 . html 
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Fig. 1. Mean field curve for B2/S2345. 



oscillators, and still life patterns. Thus unstable fixed points (Fig. 1) repre- 
sent evidence of complex behavior as mobile localizations, gliders and small 
oscillators emerging in the automaton development. Figure 2 demonstrates 
a typical evolution of B2/52345 starting with very small initial densities of 
state '1.' Mobile localizations emerge but do not survive for a long time. 



OSX Evolution 



" OSX Evolution 



OSX Evolution 



(a) 



(b) 



(c) 



Fig. 2. Snapshots of automaton configurations of 300 x 300 cells showing activity of 
particles in B2/5'2345. (a) Initial random configuration, density of state '1' is 1.2%. 
(b) Just after eight generations gliders and nucleation patterns emerge, (c) 540th 
generations. 
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A set of minimal particles, or basic periodic structures, in rule 52/5*2345 
include one glider, two oscillators (one blinker and one flip-flop configurations) , 
and one still life configuration (see Fig. 3). The still life patterns [29, 11] rep- 
resent precipitation of an abstract reaction-diffusion chemical system imitated 
by rule 52/52345 automaton. The still life blocks are not affected by their 
environment however they do affect their environment [26, 28]). Therefore the 
still life patterns can be used to build channels, or wires, for signal propaga- 
tion. 



(a) 



(b) 



■ ■ 



(c) 



(d) 



Fig. 3. Basic periodic structures in B2/S2345: (a) glider, (b) flip-flop oscillator, 
(c) blinker oscillator, and (d) still life configuration. 



2.1 Indestructible pattern in 52/52345 

Some patterns amongst still life patterns in the rule 52/52345 belong to a 
class of indestructible patterns (sometimes referred to as 'glider-proof pat- 
terns) which cannot be destroyed by any perturbation, including collisions 
with gliders. A minimal indestructible pattern, still life occupying a square of 
6x6 cells, is shown in Fig. 3d. 

The indestructible patterns discovered can be used to stop a 'super- 
nova' explosions in Life-like rules. Usually a Life-like automaton development 
started at an arbitrary configuration exhibits unlimited growth (generally re- 
lated to nucleation phenomenon [20]). 

In rule 52/52345 such an 'uncontrollable' growth can be prevented by a 
regular arrangement of indestructible patterns. Examples are shown in Fig. 4. 
In the first example (Fig. 4a) two gliders collide inside a 'box' made of still 
life patterns. The collision between the gliders lead to formation of growing 
pattern of '1' states. The propagation of the pattern is stopped by the inde- 
structible localizations (Fig. 4b) . In the second example, gliders collide outside 
the box (Fig. 4c) however interior of the box remains resting (Fig. 4c) due to 
impenetrable walls. Similarly, one can construct a colony of still life patterns 
immune to local perturbations. An example is shown in Fig. 5. 

The indestructibility exemplified above allows us to use still life patterns 
to channel information in logical circuits. 
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(a) 



(b) 



(c) 



(d) 



Fig. 4. Containment of growing pattern by indestructible localizations, (a) The 
explosion start from a reaction between two gliders, (b) final configuration of the 
container pattern, (c) Initial positions of gliders outside the box walled by indestruc- 
tible patterns, (d) Interior of the box is protected from the growing pattern. 
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(a) 



(b) 



Fig. 5. Still life colony immune to local perturbations: (a) initial configuration, 
where 'viruses' are shown are irregular patterns of '1' states, (b) final state demon- 
strates that initial local perturbations remain contained. 
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3 Computing with propagating patterns 



We built a computing scheme from channels — long areas of 'O'-state cells 
walled by still life blocks, and T-junctions 6 — sites where two or more channels 
join together. 




c 

Fig. 6. T-shaped system processing information. 



Each T-junction consists of two horizontal channels A and B (shoulders), 
acting as inputs, and a vertical channel, C, assigned as an output (Fig. 6). 
Such type of circuitry has been already used to implement XOR gate in chem- 
ical laboratory precipitating reaction-diffusion systems [3], and precipitating 
logical gates imitated in CA [26, 28]. A minimal width of each channel equals 
three widths of the still life block (Fig. 3d) and width of a glider (Fig. 3a). 



StinBBBBBBBBBBBtttttttttttttttttttttttt 
tt ■: Jt 
StinBBBBBBBBBBBtttttttttttttttttttttttt 



(a) 



(b) 




Fig. 7. Feedback channels constructed with still life patterns (fig. (a) and (c)) show 
the initial state with the empty channel and one glider respectively. The symmetric 
pattern represent value (fig. (b)), and non-symmetric pattern represent value 1 
(fig. (d)) late of glider reaction. 



Boolean values are represented by reaction of gliders, positioned initially 
in the middle of channel, value (Fig. 7 (a)), or slightly offset, value 1 (Fig. 7 
(c)). The initial positions of the gliders determine outcomes of their reaction. 



T-junction based control signals were suggested also in von Neumann [44] works. 
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Glider, corresponding to the value is transformed to a regular symmetric 
pattern, similar to frozen waves of excitation activity (Fig. 7b). Glider, rep- 
resenting signal value 1, is transformed to transversally asymmetric patterns 
(Fig. 7d). Both patterns propagate inside the channel with constant, advanc- 
ing unit of channel length per step of discrete time. 

3.1 Implementation of logic gates and beyond 

When patterns, representing values and 1, meet at T-junctions they compete 
for the output channel. Depending on initial distance between gliders, one of 
the patterns wins and propagates along the output channel. Figures 8, 9 and 11 
show final configurations of basic logical gates. 




(c) (d) 

Fig. 8. Two kinds of OR gate implementation at the Life rule B2/S'2345. Input 
binary values a and b they are represented as 'In/0' or 'In/1,' output result c is 
represented by 'Out/0' or 'Out/1.' (a) and (c) initial configurations of the gates, (d) 
and (e) final configurations of the gates. 

Figure 8 shows two implementations of OR gate. Due to different locations 
of gliders in initial configurations of gates, patterns in both implementations 
of gates are different however, results of computation are the same. Configu- 
rations of AND gate are shown in Fig. 9. 
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Fig. 9. Two kinds of AND gate implementation at the Life rule B2/S2345. (a) and 
(c) initial configurations of the gates, (d) and (e) final configurations of the gates. 



Also we can implement a delay element as shown in Fig. 10. 
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Fig. 10. Configurations of delay element for signal '0' (ab) and signal 'f ' (cd). (a) 
and (c) initial configurations, (b) and (d) final states. 
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The not gate is implemented using additional channel, where control pat- 
tern is generated, propagate and interfere with data-signal pattern. Initial and 
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final configurations of NOT gate are shown in Fig. 11. Using delay elements we 
can construct serial channels with any number of not gates (Fig. 12). Number 
of control channels growth proportionally to number of gates in the circuit. 
Of course, it could not be the most elegant and efficient way of constructing 
not gate, but useful for our purposes in i?2/5'2345. 

3.2 Majority gate 



a a i 




Fig. 13. Majority gate implementation (a A b) V (a A c) V (b A c) : (a) initial config- 
uration: majority input values In/0 (first column), and majority input values In/1 
(second column), and (b) final configurations of the majority gates. 
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Implementation of majority gate is shown in Fig. 13. The gate has three 
inputs: North, West and South channels, and one output: East channel. Three 
propagating pattern, which represent inputs, collide at the cross-junction of 
the gate. The resultant pattern is recorded at the output channel. Similarly 
gates in quantum-dot cellular automata are designed [37]. 

3.3 Implementation of binary adders 

We represent two ways of implementing partial and full binary adders. First 
we can consider design based on cascading of logic gates, then second design 
employing only not-MAJORITY gates. 



Adder via cascading of logical gates 

Conventional logic circuit and truth table of a binary half-adder are shown 
in Fig. 14a (true tables were derived from [28]). 7 Schematic representation of 
the half-adder via T-junctions between channels is shown in Fig. 14b. 



-)-^ — ' 



'•out 



a 


b 


Gout 


sum 

















1 





1 


1 








1 


1 


1 


1 






preserving 
carry oul 
or sum 



(a) 



(b) 



Fig. 14. (a) Half-adder circuit and its respective true table, and (b) scheme of 
half-adder implementation in geometrically constrained medium. 



Final configurations of patterns in a one-bit half-adder are calculated for 
all inputs values. Figure 15a shows the initial configuration of inputs and 
outputs. Resultant patterns are shown in Fig. 15b where the last patterns is 
a carry-out operation for the next half-adder. 

7 http : //uncomp . uwe . ac .uk/genaro/Dif f usion_Rule/B2-S2345678 .html 
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Fig. 15. Half adder implemented in rule B2/S'2345. Configurations represent sums 
+ 0, + 1, 1 + 0, and 1 + 1 respectively, carry-out is preserved. 
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The circuit can be extended to a full binary adder via cascading of logical 
gates [19]. Configuration of the adder, built of still life blocks, and its descrip- 
tion are shown in Fig. 16. The full adder consists of 16 T-junctions, linked 
together by channels and involve synchronization signals. 
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Fig. 16. Configuration and description of full binary adder. 



Working prototype of full adder is constructed in 1, 118 x 1, 326 cells lat- 
tice. In total, initial configuration has a population of 28,172 alive cells. The 
prototypes working cycle is 1,033 time steps with a final population of 63,662 
alive cells. A data-area of the full adder is shown in Fig. 17. 
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Fig. 17. Zoomed-in data-area of the serial implementation of full binary adder. 
Construction of adder via not-MAJORITY gates 

The second model represents a binary adder constructed of three NOT- 
MAJORITY gates and two inverters. Such type of adder appears in several pub- 
lications, particularly in construction of the arithmetical circuits in quantum- 
dot cellular automata [37]. Original version of the adder using not-majority 
gates was suggested by Minsky in his designs of artificial neural networks [31]. 

Figure 18 shows the classic circuit and the true table illustrating the dy- 
namics of this adder. This way, figure 19 represents a scheme of the adder 
to implement in 52/5*2345. The scheme highlights critical points where some 
extra gates are necessary to adjust inputs and synchronize times of collisions. 

Figure 20 presents most important stages of the full adder on 52/5*2345 
evolution space standing out delays stages and not gates. The adder is 
implemented in 1, 402x662 lattice that relates an square of 928,124 cells lattice 
with an initial population of 56,759 cells in state Final configurations of the 
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Fig. 18. Circuit and truth table of a full binary adder comprised of not-majority 
gates and inverters. 



J 




Fig. 19. Schematic diagram of a full binary adder comprised of not-majority 
gates. Delay elements are not shown. 

adder for every initial configuration of inputs are shown in Figs. 21-28 with 
a final population of 1,439 alive cells on an average of 129,923 generations. 




Fig. 20. Full binary adder designed with not-majority gates, stages and main 
circuit implementation on B2/S2345 evolution space. 




Fig. 21. Final configuration of the adder for inputs a — 0, b — and c in = 0, and 
outputs c out = and sum = 0. 




Fig. 22. Configuration of the adder for inputs a — 0, b = 1 and c in — 0, and outputs 
c ut = and sum = 1. 




Fig. 24. Configuration of the adder for inputs a = 1, b = 1 and c in — 0, and outputs 
Cout = 1 and sum = 0. 




Fig. 26. Configuration of the adder for inputs a — 0, b = 1 and a n = 1, and outputs 
Cout = 1 and sum = 0. 




Fig. 27. Configuration of the adder for inputs a = 1,6 = and c in — 1, and outputs 
Cout = 1 and sum = 0. 




Fig. 28. Configuration of the adder for inputs a = 1, b = 1 and c;„ = 1, and outputs 
Cout = 1 and sum = 1. 
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4 Conclusions 

We studied a cellular automaton model — Life-like rule i?2/S , 2345 - of a 
precipitating chemical system. We demonstrated that chaotic rule £> 2/5*2345 
supports stationary (still lifes) and mobile (gliders and propagating patterns) 
localizations. That relates another case where a chaotic rule contains non 
evident complex behavior and how such systems could have some computing 
information on its evolution [27, 32]. 

We have shown how construct basic logical gates and arithmetical circuits 
by restricting propagation of patters in a channels made of indestructible sta- 
tionary localizations. Disadvantage of the approach presented is that comput- 
ing space is geometrically constrained and not all cells of automaton lattices 
are used in computation. However, the geometrical constraining brings some 
benefits as well. Most computing circuits in Life-like automata are built using 
very complex dynamics of collisions between gliders and still lifes [9, 39, 18], 
in our case gliders are used only to 'ignite' propagating patterns in the chan- 
nels [3, 45]. In future studies we are planning to implement the computing 
architecture designed in the paper to manufacture experimental prototypes 
of precipitating chemical computers; they will be based on crystallization of 
'hot ice' [6]. 

Implementations and constructions are done in Golly system (http : / / 
golly.sourceforge.net/). Source configurations and specific initial condi- 
tion (RLE files) to reproduce the results are available in http : //uncomp . uwe . 
ac . uk/ genaro/Lif e_dc22 . html. 
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